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ABSTRACT 

We present an analysis of the globular cluster (GC) system of the nucleated dwarf elliptical galaxy 
VCC 1087 in the Virgo cluster, based on Keck/LRIS spectroscopy and archival Hubble Space Tele- 
scope/Advanced Camera for Surveys imaging. We estimate VCC 1087 hosts a total population of 
77 ± 19 GCs, which corresponds to a relatively high V-band specific frequency of 5.8 ± 1.4. The 
9 475 — z 850 color distribution of the GCs shows a blue (metal-poor) peak with a tail of redder (metal- 
rich) clusters similar in color to those seen in luminous ellipticals. The luminosity function of the 
GCs is log-normal and peaks at M™ 5 = —7.2 ± 0.3, Mj® = —8.1 ± 0.2. These peak positions are 
consistent with those found for luminous Virgo ellipticals, suggesting either the lack of, or surprisingly 
similar, dynamical destruction processes of GCs among dwarf and giant galaxies. Spectroscopy of a 
subsample of 12 GCs suggests that the GC system is old and coeval ( >10 Gyr), with a fairly broad 
metallicity distribution (-1.8 <[m/H] <-0.8). In contrast, an integrated spectrum of the underlying 
galaxy starlight reveals that its optical luminosity is dominated by metal-rich, intermediate-aged stars. 
Radial velocities of the GCs suggest rotation close to the major axis of the galaxy, and this rotation 
is dynamically significant with (v m t / <J\ os )* > 1. A compilation of the kinematics of the GC systems 
of 9 early-type galaxies shows surprising diversity in the {v mi /a\ os ) parameter for GC systems. In 
this context, the GC system of VCC 1087 exhibits the most significant rotation to velocity disper- 
sion signature. Dynamical mass modeling of the velocity dispersion profile of the GCs and galaxy 
stars suggest fairly constant mass-to-light ratios of ~ 3 out to 6.5 kpc. The present observations 
can entertain both baryonic and non-baryonic solutions, and GC velocities at larger radii would be 
most valuable with regard to this issue. Finally, we discuss the evolution of VCC 1087 in terms of 
the galaxy "harassment" scenario, and conclude that this galaxy may well be the remains of a faded, 
tidally perturbed Sc spiral. 

Subject headings: globular clusters: general - galaxies: individual (VCC 1087) - galaxies: dwarf - 
galaxies: kinematics and dynamics 



1. INTRODUCTION 

Characterized by their faint luminosities (My > —18) 
and low effective surface brightnesses (/iy c ff > 22), dwarf 
elliptical (dE) galaxies are numerically the dominant 
galaxy type in galaxy clusters (Binggeli, Sandage & Tam- 
mann 1988). dE galaxies are often subdivided into those 
with nuclei (dE,N) and those without (dE,noN), a di- 
vision with a broad luminosity dependence in the sense 
that most dEs brighter than My ~ —16 have nuclei, 
while those fainter than My ~ —12 do not (Sandage, 
Binggeli & Tamman 1985). Recently it has come to light 
that a further division may be made; one of rotating 
and non-rotating dEs (De Rijcke et al. 2001; Pedraz et 
al. 2002; Geha et al. 2002, 2003). Geha et al. (2002, 2003) 
have shown that, so far, this distinction in dE galaxy 
kinematics does not extend to their stellar populations 
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or morphologies. 

Formation scenarios for dEs must satisfy the observa- 
tional constraints that (Conselice 2004): (i) dEs pref- 
erentially exist in dense regions (even though they are 
spatially diffuse), (it) dEs exhibit both rotational and 
pressure support, and (Hi) dEs contain old, metal-poor 
stars and/or young, metal-rich stars. The presently fa- 
vored models of dE formation may be roughly distilled as 
follows: dEs are primordial objects which formed early 
on in dense environments (e.g., White & Frenk 1991) and 
possibly "squelched" by reionisation (Tully et al. 2002), 
or that dEs are objects which had a progenitor popula- 
tion (spiral or irregular galaxies) and were morphologi- 
cally transformed in the cluster environment (Mao & Mo 
1998; Moore et al. 1998; Mastropietro et al. 2005). Ob- 
servations of the velocity and spatial distributions of dEs 
in galaxy clusters suggest that they may have infallen 
'u^ui3».ar^Qatii^rrilj®04iflBiiQagi;t al. 2001). This combined 
with the fact that the stellar populations of dEs differ 
from field dwarfs (van Zee et al. 2004a), and occasionally 
show spiral structures and H I (e.g., Jerjen et al. 2000; 
Barazza et al. 2002; Graham et al. 2003; De Rijcke et 
al. 2003) supports the notion that at least some dEs were 
morphologically transformed from late- type galaxies. 
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Globular clusters (GCs) provide a different perspec- 
tive with which to investigate dE galaxy formation. 
GCs are among the oldest observable stellar populations. 
Their high binding energies makes them relatively ro- 
bust against tidal disruption, and their luminous, com- 
pact nature allows for easy separation from the galaxy 
background. These and other factors make them ideal 
tracers of the mass aggregation history of galaxies (e.g., 
West et al. 2004). Durrell et al. (1996a) obtained deep 
i?-band imaging of 1 1 Virgo dEs and found that the spe- 
cific frequencies (Sn) of their sample ranged mostly from 
3-8, similar to luminous elliptical galaxies, but consider- 
ably higher than the Sn < 1 of late-type disks. A Hub- 
ble Space Telescope (HST) /WFPC2 study by Miller et 
al. (1998) (24 Virgo, Fornax & Leo dEs) confirmed these 
findings, and found differences between the Sn of dE,Ns 
(( S N ) = 6.5 ± 1.2) and dE,noNs (( S N ) = 3.1 ± 0.5). 
Miller et al. also found that, for dE,Ns, Sn increased 
with My, but little or no such correlation was seen in the 
dE,noN class. The above results lead Miller et al. (1998) 
to argue that the progenitors of dEs were probably a het- 
erogeneous population, possibly reflected by the presence 
(or lack thereof) of dE nuclei. Lotz et al. (2004) exam- 
ined the relationship between the colors of the GC sys- 
tems and those of the galaxy envelope and nucleus using 
a larger sample of 69 Virgo, Fornax & Leo dEs imaged 
with #ST/WFPC2. These authors found that the (pho- 
tometric) metallicities of the dE GC systems correlated 
with host galaxy magnitude as Zqc oc L^ 22 , suggesting 
that the GCs "knew" about the host galaxy to which 
they would ultimately belong. It now seems clear that, 
similar to the situation for "red" GCs in luminous ellip- 
ticals, the mean colors of the "blue" GCs are correlated 
with host galaxy magnitude (Larsen et al. 2001; Strader 
et al. 2004, 2005; Peng et al. 2005). 

In this paper we present a photometric and spectro- 
scopic analysis of the GC system of VCC 1087, a nucle- 
ated dwarf elliptical (dE3,N) located in projection ^300 
kpc southwest of the cD galaxy M87 (VCC: Virgo Clus- 
ter Catalog of Binggeli, Sandage & Tamman 1985). This 
paper is organized as follows: in Section we discuss 
the observations and reductions of the data used in this 
paper. Section presents our analysis of the photomet- 
ric and spectroscopic data. In Section 0] we discuss the 
results of this study, and summarize our findings in Sec- 
tion [5] 

We note that Jerjen et al. (2004) obtain two different 
distance estimates for VCC 1087 from surface brightness 
fluctuations of (m-M) = 31.27±0.14 and 31.39±0.18, 
which correspond to My=-17.8 and -18.0 respectively. 
Unless otherwise stated, we adopted the shorter distance 
modulus which puts VCC 1087 some 2.7 Mpc beyond the 
HST Key Project distance to Virgo ((m - M) = 30.92; 
Freedman et al. 2001). 

2. THE DATA 
2.1. Imaging 

Photometry of VCC 1087 and its GC system was 
derived from Sloan-type g-band (F475W) and z-band 
(F850LP) wide-field images taken with the Advanced 
Camera for Surveys (ACS) as part of the ACS Virgo 
Cluster Survey (Cote et al. 2004). We follow the conven- 
tion of Cote et al. (2004) and refer to the ACS F475W 



and F850LP filters as 3475 and zsso respectively. 

Images were processed through the standard ACS 
pipeline. Multidrizzle was then used to combine individ- 
ual images and reject cosmic rays. Candidate GCs were 
selected as 4cr detections (matched in the 3475 and Z850 
bands) identified from 20 x 20 pixel median-subtracted 
images. Aperture photometry in a 5-pixel aperture was 
then obtained for all detected objects using the PHOT 
task in DAOPHOT II. Aperture corrections to a 10- 
pixel aperture were calculated from bright, isolated ob- 
jects in VCC 1226, VCC 1316, VCC 1978, VCC 881 and 
VCC 798 (the five brightest galaxies in the ACS/Virgo 
survey of Cote et al. 2004). These magnitudes were then 
corrected to a nominal infinite aperture using the val- 
ues of Sirianni et al. (2005). We used photometric ze- 
ropoints from Sirianni et al. to convert from the instru- 
mental magnitude scale to the AB magnitude system, 
correcting for Galactic foreground extinction using the 
reddening curves of Cardelli et al. (1989) and the DIRBE 
dust maps of Schlegel et al. (1998). Photometric com- 
pleteness limits were determined from artificial star tests. 
The 50% completeness limit in 3475 is 26.6, and in Zsso 
is 25.5. 

On the ACS images, GCs are round, compact objects; 
they can be separated from the majority of foreground 
stars and background galaxies using shape, color and 
size information. GC candidates were selected as ob- 
jects detected in both 17475 and Zs50 images to within 0.1 
arcsec tolerance, with roundness and sharpness parame- 
ters (from DAOPHOT II) in the intervals -0.5< round- 
ness <0.5 and 0.55< sharpness <0.90 respectively. A 
color range of 0.5 < (3475 — z 850 ) < 2.0 was adopted 
for GCs, which encompasses age and metallicity combi- 
nations of 1 Gyr, [m/H]=-2.25 to 15 Gyr, [m/H]=0.67 
(Maraston 1998). Objects with 3475 < 19 and z 850 < 18 
were deemed too luminous to be true GCs. These 
bright limits are roughly two magnitudes brighter than 
the most luminous Milky Way GC (u> Cen) at the dis- 
tance of VCC 1087. Sizes for the GC candidates were 
measured using ISHAPE (Larsen 1999), and this infor- 
mation was used to expunge objects outside the range 
of half-light radius (nj expected for Milky Way GCs 
(0.5 < Hi(pc) < 12). Finally, the GC sample was visu- 
ally inspected for any obvious non-GC properties (e.g., 
visible structures, blends) and these were culled from the 
final sample. 

For g 475 < 24, GCs are easily distinguished from galax- 
ies since they appear more compact than galaxies for a 
given magnitude. Indeed, from our spectroscopy (see 
SectionE2J), all the GCs identified on the HST/ ACS im- 
ages were found to be bona fide clusters (the faintest 
being 3475 = 23.7). Therefore, we conclude that the 
contamination rate for cluster candidates brighter than 
this (approximately the turn-over (TO) in the GCLF) 
is consistent with zero. To quantify the contamination 
from background galaxies at magnitudes fainter than 
(7475 = 24, a master background file was constructed 
from the ACS/Virgo survey dwarf galaxy fields at large 
radii as discussed in Strader et al. (2005). From a region 
totalling 41 sqr. arcmin, we identify 37 objects which 
pass our criteria for GCs where we expect none (although 
note that, for the most luminous dEs such as VCC 1087, 
the GC systems do extend to the edge of the ACS field as 
confirmed by our spectroscopy, e.g., Figure Scaling 
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to one ACS pointing, we find that ~ 10 objects may be 
background galaxies. We detect 68 GCs associated with 
VCC 1087. Based upon the incompleteness at the faint 
end of the GC luminosity function (GCLF) (Section 0, 
and correcting for incomplete areal coverage we estimate 
we have identified ~ 88% of the total GC population in 
this galaxy. Therefore we estimate a total GC popula- 
tion of 77 GCs. Our estimate of the contamination due 
to background galaxies suggests that we may be over- 
estimating the GC population by ~ 13%. Rather than 
correcting for the total number of GCs, we fold this ab- 
solute uncertainty into our population estimate to yield 
77 ± 19 GCs. This number excludes the nucleus of the 
galaxy (which meets all of our GC selection criteria.) 
Adopting My =-17.8 for VCC 1087 (Jerjen et al. 2004), 
we obtain a specific frequency of 5.8 ± 1.4. 

2.2. Spectroscopy 

Optical spectra were obtained for 12 GC candidates 
iden tified from the 17475 and zsso HST/ ACS images (Sec- 
tion I2.1|l during the nights of 24-25 April, 2004, us- 
ing the Low Resolution Imaging Spectrograph (LRIS) 
(Oke et al. 1995) on the Keck I telescope. A total of 
15 1800s exposures were obtained (totaling 7.5 hours) 
through a 0.8 arcsec slitmask. We also obtained a 
300s, 0.8 arcsec long-slit spectrum of VCC 1087 ori- 
ented along the major axis of the galaxy (PA. =104°). 
A flux standard (BD+33°2642) taken from Oke (1990) 
was observed with a longslit to allow first-order contin- 
uum corrections, along with a number of radial velocity 
and Lick/IDS standard stars (from Worthey et al. 1994). 
Bias frames, dome/sky flats and Cadmium- Neon- Argon 
arcs were taken at various points during the run for cal- 
ibration purposes. The seeing varied from 0.6-0.9 arcsec 
during the two-night run. 

We obtained simultaneous blue and red spectra 
through the use of a dichroic which split the beam at 
6800A. On the blue side, a 600 lines/mm grism blazed at 
4000A was used to yield an effective wavelength range of 
3300-5900 A and a resolution of ~3A (FWHM). This ef- 
fective wavelength range varied slightly between objects 
due to the varying slitlet positions on the slitmask. On 
the red side, a 831 lines/mm grating was used, blazed at 
8200A to cover the Ca II triplet region. The analysis of 
the Ca II triplet data will be discussed elsewhere. 

The reduction of these data was performed with 
IRAF 3 . All science images were bias-subtracted and then 
flat-fielded with twilight sky flats which were co-added 
and normalized. Spectra were optimally extracted (see 
Home 1986) and wavelength calibrated with solutions 
obtained from the arc exposures. Wavelength residuals 
of 0.2 A were typical. The zeropoints of the wavelength 
calibration were checked against skylines in the back- 
ground spectra, and were in some cases adjusted by a 
few tenths of an Angstrom. These wavelength solutions 
were further adjusted as described below. The spectra 
were then combined with cosmic-ray rejection, and flux 
calibrated using the response function derived from our 
flux standard. Error spectra were extracted in parallel 

3 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



with the object spectra to facilitate accurate estimation 
of uncertainties in our analysis. 

The velocity dispersion of the VCC 1087 GC system 
was anticipated to be similar to that of the galaxy itself 
at ~ 40 km s _1 (Geha et al. 2003). This required that 
particular care was taken to characterize the uncertain- 
ties in our radial velocity measurements in this study. 
The theoretical velocity accuracy for our spectral reso- 
lution is ~ 20 km s _1 at 4500A. Radial velocities for 
the GCs were derived using the MOVEL algorithm de- 
scribed in Gonzalez (1993). This is an improvement over 
the classic Fourier quotient method (Sargent et al. 1977) 
which determines radial velocities and velocity disper- 
sions simultaneously Although our spectral resolution 
did not allow to us measure velocity dispersions for the 
GCs themselves, we used this approach since experimen- 
tation showed that it provides more robust velocity deter- 
minations than classical cross-correlation procedures. As 
reference template spectra for this method, we employed 
a number of Lick stars observed during the run. To 
prevent overestimating the uncertainties in the method, 
which can appear when velocity dispersion solutions are 
close to zero, the GC spectra were broadened by con- 
volving them with a Gaussian of a = 50 km s _1 . Uncer- 
tainties in the radial velocities were computed from 100 
Montecarlo simulations using the error spectra. 

After measuring velocities, small deviations from 
the linear wavelength scale were corrected by using 
the MILES spectral energy distributions (Vazdekis et 
al. 2005, in preparation) as rest-frame wavelength cali- 
brators. Each GC spectrum was cross-correlated with an 
appropriate MILES template in ~100 small, overlapping 
spectral regions in the dispersion direction. The derived 
offsets (around zero but not necessarily null) were fitted 
as a function of wavelength, thus providing a low-order, 
polynomial correction to the prior wavelength calibration 
that was applied to the deredshifted data. In this way, all 
the final GC spectra resulted in having an homogeneous, 
rest-frame, linear wavelength scale. These velocities are 
listed in Tabled From Tabled the mean velocity uncer- 
tainty is 22 km s , in good agreement with the expected 
value. Example spectra are shown in Figure ^ 

Lick indices were measured from the spectra using the 
definitions given in Trager et al. (2000) and Worthey & 
Ottaviani (1997). Index uncertainties were determined 
from the error spectra. Small additive offsets where ap- 
plied to the measured indices obtained from our Lick 
standard star observations, and were in general smaller 
than the index uncertainties. 
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Fig. 1. — Comparison of two VCC 1087 GC spectra one 
of intermediate metallicity (GC#12; upper), and one metal- 
poor (GC#5; lower). Magnesium hydride (~5170A) and the 
G-band (~4300A) are stronger in the more metal-rich GC. 
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Fig. 2.— HST/ACS g 4 rs image of VCC 1087. GC candi- 
dates identified from our HST photometry are indicated by 
small circles, those for which we have spectra and are posi- 
tively identified as GCs are double-circled. The bright nucleus 
of the galaxy can be seen in the center of VCC 1087, adjacent 
to the most central GC. 



3. ANALYSIS 

3.1. Photometric Properties of the GC System 

The spatial distribution of the GC candidates on the 
5475 image is shown in Figure [21 The spectroscopically 
confirmed GCs are also indicated (Table^J . The GCs are 
concentrated towards the center of the galaxy, and ap- 
pear spatially elongated in the NE-SW sense. This gen- 
eral impression is supported by the azimuthally-averaged 
distribution of the GCs (Figure[3J). A sinusoidal fit to the 
unbinned data gives a position angle (P. A., measured E 
through N) of 85 ± 49°. Unfortunately the small sample 
size does not allow for stronger constraints on the major 
axis P. A. of the GC system, but the P.A.'s of the GCs 
and the galaxy isophotes (104°) are consistent with each 
other. The major- to minor-axis ratio of the GC spatial 
distribution gives an approximate ellipticity of e ~ 0.3, 
which is similar the galaxy ellipticity (e = 0.26) beyond 
1 kpc. 
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Fig. 3. — Angular distributions of GCs associated with 
VCC 1087 interior to 1.4 arcmin (the radius from the galaxy 
center which permits full azimuthal coverage; shaded his- 
togram) and for the full sample of GCs (dashed line). The 
vertical dashed line represents the position angle of the major 
axis of VCC 1087 (P.A -104°), the short dashed line indicates 
the position of the minor axis. 

The nucleus of the galaxy is well resolved, and we mea- 
sure 3475 = 20.22 ± 0.01, 3475 - ^sso = 1-30 ± 0.01 (Ta- 
bled . We have also measured the envelope colors of the 
galaxy, and obtain 3475 — Zsso = 1-38 ± 0.01, indicat- 
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ing that the envelope of the galaxy is marginally redder 
than its nucleus 4 . This is in agreement with the HST 
studies of De Propris et al. (2005) and Lotz et al. (2004), 
but is in the opposite sense to the conclusions of ground- 
based work (e.g., Rakos & Schombert 2004; Caldwell & 
Bothun 1987). Presumably this disagreement stems from 
the fact that the ground-based work failed to resolve the 
dE nuclei. 




Fig. 4. — Color-magnitude diagram for all GC candidates 
from HST/ ACS images (open circles) and those for which we 
have obtained spectra (filled circles). The horizontal dashed 
line (upper panel) shows the 50% completeness limits of these 
data. 

The color-magnitude distribution for the confirmed 
and candidate GCs is shown in Figure The GCs show 
a blue peak (metal-poor clusters, henceforth MPCs), 
and a possible tail of red (metal-rich clusters, henceforth 
MRCs) GCs. A straight mean of the color distribution 
gives 3475 — ^850 = 0-98 with a dispersion of 0.22 mag. 
Peak locations for the MPCs and (putative) MRCs were 
determined using the Bayesian code Nmix (Richardson 
& Green 1997), which implements mixture modeling of 
heteroscedastic normal distributions with the number of 
subpopulations a free parameter. For VCC 1087, the pos- 
terior distribution of this parameter peaked at two, sug- 
gesting two groups in the color distribution with peaks 
at 3475 - ^850 = 0.90 and 1.17 for the MPCs and MRCs 
respectively. 

The red objects (3475 — 0g5O ~ 1-5) near to our 50% 
completeness limit are possibly background galaxies (this 
statement can only be confirmed spectroscopically). To 
check whether these objects are driving this bimodal sig- 
nature, we cut the CMD at 3475 = 25 and re-ran the 
code. Again, two populations were preferred with peaks 
at 3475 - z 850 = 0.90 and 1.16. Similarly, a KMM 
test (Ashman, Bird & Zepf 1994) returned two peaks 
a t 3475 — z 850 — 0-89 and 1.16, although with a fairly 



low significance (p = 1.8). These peak positions are con- 
sistent with the galaxy luminosity- GC color relations 
derived by Strader et al. (2005). We tentatively identify 
the objects with 3475 ~ 24 and 3475 — zsso ~ 1-1 as 
MRCs, three of which, confirmed by our spectroscopy, 
are consistent with this interpretation (e.g., Figure EJ- 
Note that, in the subsequent analysis, we have not sep- 
arated out the blue and red subpopulations due to the 
small number of clusters. 

The radial density profile of the GC system is com- 
pared to the galaxy surface brightness profile in Fig- 
ure El We fitted the GC profile with the functional form 
p (x r~ a . Excluding the innermost bin (see below), we 
find a = — 1.1 ± 0.4, where the uncertainty comes from 
a bootstrap of the density profile. A power-law index 
of a = — 1.1 ± 0.4 is relatively shallow, more charac- 
teristic of the GC systems of cD galaxies than those of 
low-luminosity ellipticals (e.g., Harris et al. 2004). Simi- 
larly shallow density profiles have been suggested for the 
GC systems of dEs previously (Durrell et al 1996a; Lotz 
et al. 2001). The galaxy starlight is well fit by a Sersic 
(1968) profile with n = 1.4 (Geha et al. 2003), where 
n = 1 corresponds to an exponential profile, and n = 4 
corresponds to an r 1 / 4 law. 

There is the suggestion of a "core" in the GC density 
profile near the galaxy effective radius, where there is a 
deficit of clusters compared to the power-law slope. Lotz 
et al. (2001) discussed the existence of such cores in the 
context of the formation of dE nuclei through the dy- 
namical destruction of GCs at small radii. Could these 
"missing" GCs be the cause of the break in the den- 
sity profile? For GCs at the turn-over of the GCLF (see 
below), Lnucieus/icc ~ 40, i.e., 40 typical GCs would 
need to coalesce to give the observed nuclear luminos- 
ity. These additional GCs (added to the innermost ra- 
dial bin) are indicated by the open square in Figure [S] 
(for simplicity, we assume that these "missing" GCs oc- 
cupied the same radial distribution as those present in 
the inner bin.) The position of these GCs lie on an ex- 
trapolation of the power-law fit in Figure [3] This would 
seem to suggest that the existence of a core may indicate 
dynamical destruction processes. 

However, since the timescale for dynamical friction is 
expected to scale as M^, one would expect that the 
most massive GCs decay first (for a given spatial density 
profile) . Drawing clusters from the bright end of the ob- 
served GCLF (FigureEJ), ~ 10 GCs are required to equal 
the luminosity of the galaxy nucleus (the number of lu- 
minous GCs required is clearly dependent upon the exact 
shape of the "zero-age" luminosity function, but presum- 
ably the bright end has not become brighter with time). 
This is indicated by the open circle in Figure whose 
position is ~ 4 a away from the extrapolated profile fit. 
From these very simplistic arguments, unless somehow 
turn-over mass clusters are preferentially disrupted, the 
luminosity of the nucleus and the core in the inner parts 
of GC surface density of VCC 1087 cannot both be ex- 
plained purely due to the action of dynamical friction. 



4 It was unnecessary to mask out the nucleus in measuring the 
envelope colors since it contributes only ~ 1% to the total inte- 
grated light of the galaxy. 
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Fig. 5.— Radial density profile of the VCC 1087 GC system 
(filled circles with Poisson error bars) compared to the sur- 
face brightness profile of the galaxy (arbitrarily normalized) . 
The short dashed lines indicates a power-law fit to the GC 
data with exponent a = —1.1, excluding the innermost bin. 
The vertical dotted line indicates the effective radius of the 
galaxy at ~1.2 kpc. The open square shows the expected sur- 
face density of GCs "corrected" for the effects of dynamical 
friction for GCs at the turnover of the GCLF. The open circle 
shows the same but for GCs at the bright end of the GCLF 
(see text). 




Fig. 6. — Luminosity functions for the VCC 1087 GC system 
in g 475 (top panel) and z sso ■ Curves are maximum-likelihood 
estimates for the peaks and dispersions of the GCLFs assum- 
ing Gaussian distributions (short-dashed curves : the peak 
and dispersion are free parameters; long-dashed curves : the 
dispersions are fixed at cr S475 = 1.1 and cr Z 850 = 1.0). The 
vertical dashed lines represent the 50% completeness limit of 
the photometry. 

The luminosity functions for the VCC 1087 GCs in 
the g 475 and z §50 bands are shown in Figure [S] The 
50% completeness limits of the photometry suggest that 
we detect GCs down to the faint limit of the GCLF. 
We fitted both Gaussian and ts functions to the un- 
binned magnitude data using a slightly modified ver- 
sion of the maximum-likelihood code of Seeker & Har- 
ris (1993). Within the uncertainties of the fits, both 
functions yield indistinguishable means and dispersions 
(assuming a t 5 — 0.775cr ff ), and we give the more fa- 
miliar Gaussian best-fit parameters. Leaving the Gaus- 
sian means and dispersions as free parameters in the 
fit, we obtain g\% = 24.1 ± 0.3, cr g475 = 1.4 ± 0.3 and 
= 23.2 ± 0.2, a Z850 = 1.26 ±0.3. 

These fits are shown in Figure H3 Inspection of the 
figure suggests that the dispersions of the GCLF, par- 
ticularly in the case of 3475, may be overestimated. Ex- 
perimentation with the Seeker & Harris (1993) code sug- 
gested that for small N, a can become artificially large. 
The straight standard deviations of the (7475 and z^5o 
GCLFs (not accounting for incompleteness) are 1.1 and 
1.0 respectively, and we fixed these in the code to assess 
the impact on the derived GCLF peak positions. g™ 5 be- 
comes 0.06 mag brighter, while z^ becomes 0.15 mag 
brighter. Both these shifts are within the uncertainties, 
and therefore we adopt the original maximum-likelihood 
parameters. 

At our adopted distance to VCC 1087, the appar- 
ent turnover magnitudes we derive correspond to abso- 
lute magnitudes of M™ 5 = -7.2 ± 0.3 and M™ = 
—8.1 ± 0.2. These turn-over magnitudes are consistent 
with that found by Strader et al. (2005), who constructed 
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a GCLF from 37 individual Virgo dEs. These authors 
found Mj° = -7.21 ± 0.14 and M 2 TO = -8.14 ± 0.15, 

y 475 Z 850 

in which our photometry for VCC 1087 was included. 

Interestingly, from the mean GCLF turnover mag- 
nitudes of three giant Virgo ellipticals (NGC 4486, 
NGC 4472 & NGC 4649), Strader et al. (2005) found 
Mj° 5 = -7.07 ± 0.15 and M™ = -8.19 ± 0.15, 
i.e., the Virgo dEs (including VCC 1087) have turnover- 
magnitudes consistent with luminous Virgo ellipticals. 
This consistency, particularly in z g 50 for which color dif- 
ferences in the systems are expected to have little effect 
on the GCLF peak positions, seems to suggest that ei- 
ther dynamical destruction of GCs in dEs and luminous 
Es operates very similarly, or that it is not a dominant 
process in shaping the mass function of GC systems. Dy- 
namical friction (Chandrasekhar 1943) and disk/bulge 
shocking (Fall & Rees 1977) are expected to operate on 
the high- and low-mass ends of the GC mass function 
respectively. These dynamical destruction processes are 
expected to be functions of the dynamical structure and 
gravitational potential of the host galaxy (for a given GC 
mass function and spatial distribution), which makes it 
unlikely that dwarf and giant galaxies will destroy their 
GCs similarly. Processes which may act against dynam- 
ical friction are discussed in Lotz et al. (2001), and in- 
clude high dark matter halo densities and/or larger halo 
core radii with respect to the baryons, and kinetic heat- 
ing mechanisms such as that which may come from GC 
scattering from black holes or through tidal stirring. The 
affects of disk/bulge shocking could be ameliorated by 
placing the GCs on primarily tangential orbits, thereby 
avoiding the galaxy centers. These possibilities would be 
a fruitful area of exploration through numerical simula- 
tion. 



3.2. Kinematics 

We obtain a biweight location (Beers et al. 1990) veloc- 
ity for the GC system of 686 ±24 km s _1 . This compares 
well with the canonical NED value (Simien & Prugnicl 
2002) of 675 ± 12 km s , and with the velocity derived 
from our integrated spectrum of the galaxy (646 ± 30 km 
s _1 ). We have looked for rotation in the GC system by 
performing linear weighted fits to the GC velocities as 
a function of position angle. The GCs show a maximal 
rotation signature about P.A.=127±18°, with 87 ± 29 
km s _1 arcmin -1 . Over the radial extent of these data 
(i.e., measured from the center of the galaxy to the out- 
ermost GC) this corresponds to V«^ = 104 ± 35 km s _1 , 
assuming solid-body rotation. The linear fit to the GC 
velocities about P.A.= 127° is shown in Figure [7| The 
P.A. of the rotation axis of the GCs appears to be close 
to parallel to the P.A. of the galaxy isophotes (104°) and 
the P.A. of the GCs (85 ± 49°). This alignment of the 
rotation axis is in a similar sense to that found by Puzia 
et al. (2000) for the GC system of the dE NGC 3115 
DW1. The GCs in VCC 1087 appear to rotate along an 
axis of 37 ± 18°, which is close to the minor axis P.A. of 
VCC 1087's isophotes (~ 14°). 




0.5 

X [arcmin] 

Fig. 7. — Radial velocities of the GCs (corrected for a sys- 
temic velocity of 686 km s _1 ) as a function of projected dis- 
tance along a P.A. of 127°. The dashed line represents a 
weighted linear fit to these data, suggesting rotation in the 
GC system of 52±18 km s _1 . The dotted line indicates an 
unweighted linear fit to these data and yields comparable ro- 
tation. 

The biweight scale of the velocities (velocity disper- 
sion, cr los ) of the VCC 1087 GCs is 49 ± 15 km s" 1 . To 
obtain a "true" velocity dispersion for the GCs (denoted 
<j[ os ) we corrected for our intrinsic velocity uncertainties 
(subtracting in quadrature) and for rotation in the GC 
system. This yields cr[ — 29 ± 6 km s _1 at a mean 
galactocentric radius of 35 arcsec (~ 3 kpc). Similar to 
such considerations for elliptical galaxies, one may ask 
whether the rotation in the GC system is dynamically 
significant by considering an isotropic oblate rotator with 
ellipticity e (e.g., Binney 1978): 



{y<j>i & los) observed fa^fo" los)observed 



model 



£ /(l-e)]V2 



(1) 



For a system with e = 0.3 (approximately the elliptic- 
ity of both the galaxy the GC system), {v^/oios) model = 

0.65. For the VCC 1087 GCs we find: (z^/o^Jobscrvcd = 
3.6±1.8, giving {v^/ <j\ os )* = 5.5±2.6. This result is sug- 
gestive that the GC system is supported by rotation. 

We have investigated where the VCC 1087 GC system 
lies with respect to other GC systems in the ellipticity— 
(v^/aios) plane. There is surprisingly little information 
on the azimuthal distributions GC systems in the litera- 
ture, precluding an accurate knowledge of e. We have as- 
sumed that the ellipticity of the host galaxy reflects that 
of the GC system (and, where they have been separated, 
that of the both the MPC and MRC subpopulations). 
There is some observational support that this assumption 
holds for the MRCs (Dirsch et al. 2005; Lee et al. 1998) 
although for the MPCs this is questionable (Kundu & 
Whitmore 1998, but see Section EHJ . In the Milky Way 
GC system, the (ensemble) halo GCs are broadly spher- 
ical, and the metal-rich GCs show flattening about the 
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Galactic plane (e.g., Borkova & Marsakov 2000). These 
considerations notwithstanding, we plot the (v^/ai os ) ra- 
tios for 9 GC systems in Figure |S| using the data and 
sources given in Tabled 

Only the GC systems for early-type galaxies are shown 
in Figure 03 since i^/cios is less meaningful for disk sys- 
tems (note that we classify NGC 4594-the Sombrero-as 
an SO here.) The GC systems of the four E galaxies lie 
on or below the curve for rotationally flattened systems 
(i.e., velocity dispersion dominates rotation), whereas the 
positions of the two dE galaxies suggests that rotation 
becomes increasingly important for these systems. The 
kinematics of the GC systems which have been separated 
into MPC and MRC subpopulations show significant di- 
versity, with neither subpopulation biased to either rota- 
tional or pressure support. In all cases, if the MPCs are 
spherically distributed, similar to the case for the Milky 
Way halo clusters, the relative importance of rotation 
would increase. 

It is interesting to compare the kinematics of the 
GC systems of NGC 5128 and NGC 4472 (at e - 0.2 
in Figure [SJ, both of which straddle the (v^ / <7\ os )* = 
0.5 curve suggestive of velocity anisotropy. Consider- 
ing the differences in luminosity and environment be- 
tween these two systems, the similarities in their kine- 
matics are quite surprising (L, V^a^/L, Vnsi28 ~ 2; 
NGC 4472 resides in the Virgo cluster, avirgo ~ 700 
km s -1 ; Binggeli et al. (1993), NGC 5128 dominates a 
small group, fTN5i28group ~ 120 km s _1 ; van den Bergh 
(2000)). Whatever the formation histories of these dis- 
parate galaxies may be, they have clearly achieved quite 
similar global kinematic signatures in their GC systems. 
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Fig. 8. — The ratio of rotation to velocity dispersion 
(i^/cios) °f the GC systems of spheroids versus host galaxy 
ellipticity. Connected filled and open circles represent MPC 
and MRC subpopulations, single filled symbols represent 
systems in which subpopulations (if any) have not been 
separated. Circles indicate elliptical galaxies, squares-SOs, 
triangles-dEs. The arrow shows NGC 3115 DW1 after cor- 
recting its GC velocity dispersion for observational errors of 
50 km s _1 . The solid curve indicates t^/cios (Equation 
for an oblate rotator of ellipticity e flattened by rotation 
(i.e., (v^/aios)* = 1). Curves for («^/cr los )* =0.3, 0.5, 2 & 
3 are also shown (long-dashed, dotted, short-dashed and dot- 
dashed curves respectively) . 



3.3. Mass-to-Light Ratios 

GCs represent one of the few probes of the underlying 
mass distribution in gas-poor galaxies beyond the stellar 
light. In fact, they are particularly useful in the case of 
dE galaxies which generally have no luminous satellites 
and often have low-surface brightness envelopes which 
limit the utility of studies of their integrated starlight. 

In order to model the radial mass (and mass-to-light 
ratio, Ty) profile of VCC 1087, we first made the as- 
sumption that the VCC 1087 GCs may be regarded as 
an extension of the galaxy starlight. That is to say, we 
explicitly assume that the stellar envelope and GCs have 
the same orbital and spatial distributions. In so doing, 
we fit dynamical models to the velocity dispersion pro- 
files of the galaxy stars (from Geha et al. 2003) as a 
function of radius by solving the spherically symmetric 
Jeans equation (see e.g., Geha et al. 2002; Gebhardt & 
Fischer 1995; van der Marel 2004): 

GM(r) 



where p is the mass density derived from the galaxy sur- 
face brightness profile and o~ r is the radial velocity disper- 
sion. The parameter describes the velocity dispersion 
anisotropy, which we assume is constant with radius. 
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of the galaxy. Assuming an isothermal-like potential, C 
is given by (Evans et al. 2003): 
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Fig. 9. — Velocity dispersion of the stellar light (solid circles) 
and GCs (open square) of VCC 1087 as a function of galacto- 
centric radius. Short and long dashed lines indicate isotropic 
and anisotropic models of constant Ty fit to the stellar data 
only (see text). A constant of 10 arcsec has been added to 
the a;-axis for display purposes such that log (R+l kpc)=3 
corresponds to the center of VCC 1087. 

The velocity dispersion of the galaxy stars and GCs as 
a function of galactocentric radius is shown in Figure 
In the figure, two models of constant Ty are shown, fit 
to the stellar data only. One is for an isotropic velocity 
dispersion ((3 = 0) and one is radially anisotropic (j3 = 
0.4). A \ 2 test indicates that the anisotropic model is 
slightly preferred over an isotropic velocity dispersion. In 
either case, Figure suggests Ty ~ 3 as traced by the 
GCs at a mean radius of ~ 3 kpc. 

The right-hand side of Equation [21 only describes the 
contribution of "pressure" to the total mass of the sys- 
tem. This is appropriate for the galaxy stars which show 
no sign of rotation (Geha et al. 2003), but is incorrect for 
the GCs for which rotation appears dynamically signifi- 
cant (Section 13. 2|) . Moreover, with only 12 clusters dis- 
tributed over a wide range of radii, the validity of using 
the Jeans equation at all must be questioned. Perhaps 
more appropriately, we can use statistical estimators to 
estimate the mass traced by the GCs which will also al- 
low for a check on the above results. 

Evans et al. (2003) introduced the tracer mass esti- 
mator (TME), which is a generalization of the projected 
mass estimator for "test particles" (Heisler et al. 1985). 
We have used the TME in order to estimate the total 
gravitating mass enclosed within concentric shells out to 
the projected radius of the outermost GC. The TME 
takes the form: 

M(R < i? max ) = — v l R i w * ( 3 ) 

i 

where N is the number of test particles (in this case 
GCs), Vi is the line-of-sight velocity of the ith GC and 
Ri is the projected distance of the zth GC from the center 
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here, r m and r out are the inner and outer radii at which 
the population density falls off as r -7 . Note that we 
have included a weighting term, ti^, in Equation @ to 
prevent undue influence being given to GCs with large 
velocity uncertainties, 



J2m = J2 a 4 



N 



(5) 



where a v> i is the velocity uncertainty of the ith GC. 

We construct three radial bins for which (r ln , r out ) are 
(0,2), (2,4) & (4,9) kpc respectively. We also set 7 = 2.1 
for all three bins. The density profile of the GC appears 
to flatten in the innermost bin (Figure |SJ ; reducing 7 
by 50% leads to - 65% reduction in Ty. The "light" 
in Ty comes from the V-band surface brightness profile 
of Geha et al. (2003) (Figure [5}, scaled to integrate to 
My=-17.8 (Jerjen et al. 2004; Barazza et al. 2003). 

Values for Ty estimated using the TME are shown as 
solid circles in Figure ITU1 Using this approach, we find 
Ty ~ 3 which is consistent with that estimated from 
Equation [21 However, as noted previously, the GC sys- 
tem shows significant rotation. The curve in Figure 1101 
indicates Ty(r), assuming solid-body rotation in the GC 
system, with M(r) = v^(r)r/G. As advocated by Evans 

et al. (2003), M pre ssure+Mratation=Mtotai, and our corre- 
sponding values for Ty are shown by the open circles in 
Figure El 

The principle uncertainties in the mass estimates stem 
from the unknown orbital configuration of the GCs, the 
systemic velocity of the GC system, and the number of 
clusters in the sample. The uncertainty in the distance 
to VCC 1087 corresponds to errors in the mass estimates 
of only ~ 5%. Within these relatively large uncertain- 
ties, even with the inclusion of Mrotation, Ty is flat with 
radius, and is consistent with that estimated from the 
stellar light at small radii. For a solar metallicity, 5 Gyr 
old stellar population (Section 00} the Bruzual & Char- 
lot (2003) models predict Ty of 3 (2) for a Salpeter 
(Chabrier) initial mass function (IMF). With the as- 
sumption of an entirely old (12 Gyr), solar metallicity 
stellar population, the Bruzual & Chariot (2003) predict 
Ty=7 (4) for a Salpeter (Chabrier) IMF. Similar val- 
ues are found for the Maraston (1998) stellar population 
models. Therefore, Ty measured from the GCs out to 
~ 6.5 kpc is consistent with that expected for a purely 
baryonic, old stellar population. Kinematics for a larger 
number of GCs at large radii would be extremely valu- 
able for further investigation of the mass distribution of 
VCC 1087. 
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Fig. 10. — Mass-to-light ratio as a function of radius in 
VCC 1087. Solid circles are the predictions assuming purely 
pressure support in the tracer mass estimator (Evans et 
al. 2003), open circles include the addition of mass from the 
rotation in the GC system (dashed curve). Filled and open 
diamonds show the predictions for the stellar light (Geha et 
al. 2003) using Equation |5] and assuming an anisotropy pa- 
rameter (/3) of 0.4 and respectively (see text). 



3.4. Ages and Abundances 

We have investigated the stellar populations of the 
GC system by using Lick indices (Table [3J in conjunc- 
tion with the simple stellar population (SSP) models 
of Thomas, Maraston & Bender (2003), and Thomas, 
Maraston & Korn (2004). These models (hereafter col- 
lectively referred to as TMK04), attempt to account 
for non-solar abundance ratios in simple stellar popu- 
lations by using synthetic atmosphere models to pre- 
dict Lick index "response functions" for varying abun- 
dance ratios. Note that the assumption of an SSP prob- 
ably breaks down for a dE as luminous as VCC 1087 
(Ly ~ 1.2 x 10 9 L Q ). Therefore, the following derived 
stellar population parameters of VCC 1087 itself are 
strictly luminosity-weighted properties only. 

In Figure im we show an example index-index diagnos- 
tic plot of the Lick indices of the GCs compared to the 
TMK04 models. We have chosen the [a/Fe]=0.3 mod- 
els here (see below). Qualitatively, the GCs generally 
appear old, and have metallicities less than solar. The 
galaxy starlight, on the other hand, is more metal-rich, 
and adopting the [a/Fe]=0.0 models for this galaxy, lies 
between the 4 and 5 Gyr model lines. Also shown in 
FigurelTTI are Lick indices for 40 Galactic GCs measured 
from the library of Schiavon et al. (2005). These GCs 
form a tight sequence, reflecting the high quality of these 
data. However, they are clearly offset from the models 
(as are the VCC 1087 data, in the mean), an offset which 
is significantly larger than the uncertainties in correcting 
to the Lick system. Investigating the origins of this offset 
(seen in the other Balmer indices to lesser or greater de- 
grees) is beyond the scope of this paper, but clearly needs 



to be resolved before confidence can be put into the abso- 
lute age-scale of SSP models (e.g., Lee & Worthey 2005). 
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Fig. 11.— The [MgFe]-H 7P indices of the VCC 1087 GCs 
(solid circles) compared to the stellar population models of 
TMK04. Lick indices measured from the Galactic GC sample 
of Schiavon et al. (2005) are shown as open circles. The heavy 
line indicates an age of 15 Gyr, the dashed line indicates an 
age of 3 Gyr. A line of constant solar metallicity is also shown. 
The GCs are appear generally old within the uncertainties, 
VCC 1087 itself (open triangle) appears more metal-rich and 
younger. 

To place these impressions on a quantitative foot- 
ing, metallicities (here denoted [m/H]), ages and a- 
abundances (here denoted [E/Fe] - see Proctor & San- 
som 2002) were derived for the GCs and galaxy starlight 
from their Lick indices by performing multivariate fits to 
the TMK04 models. This technique has been described 
extensively in Proctor & Sansom (2002) and Proctor 
et al. (2004). In general, 10-15 indices were fit simul- 
taneously, and we show the resulting best-fit ages and 
metallicities using this procedure in FigurelT2l The same 
procedure was carried out for the Schiavon et al. (2005) 
data. The VCC 1087 GCs exhibit a range of metallici- 
ties (—1.8 <[m/H] < — 0.8) with the three most metal- 
rich GCs entering into the domain of the red peak of 
the Galactic GC system. Their position in Figure IT2l is 
consistent with that shown in Figure i.e., at or near 
the separation of the two subpopulations. Note that the 
Schiavon et al. library comprises of ~ 25% of the known 
GCs in the Milky Way, whereas our spectroscopic sample 
covers the brightest ~ 15% of the VCC 1087 GC system. 

Figure suggests that the VCC 1087 GCs are all old, 
comparable to the Galactic GCs ages, and within the un- 
certainties (typically 3 Gyr) are coeval. By way of con- 
trast, the luminosity-weighted properties of VCC 1087 
itself are quite different. It is metal-rich ([m/H]= —0.1), 
lying between the metal-rich Galactic GCs NGC 6528 
and NGC 6553 in Figure lTTI and intermediate-aged (~ 4 
Gyr). This is consistent with the findings of Geha et 
al. (2003) who also used Lick indices and the TMK04 



Virgo Dwarf Globular Clusters 



11 



models. 




5 10 
No. of GCs 



Fig. 12.— Ages and metallicities of the VCC 1087 GCs (solid 
circles) predicted using the TMK04 stellar population models. 
Also shown is the luminosity-weighted age and metallicity of 
the integrated starlight of VCC 1087 (open square) , and of a 
sample of 40 Galactic GCs from Schiavon et al. (2005) (open 
circles). The VCC 1087 clusters are old and coeval, with a 
range of metallicities. The shaded and dashed histograms are 
for the VCC 1087 and Galactic samples respectively. 




Fig. 13. — Predicted a-element abundances of the VCC 1087 
GCs using the TMK04 models (solid circles) and the galaxy 
starlight (open square) versus metallicity. The Galactic 
GCs of Schiavon et al. (2005) are also shown (open circles) 
The dashed line (dotted lines are the rms scatter) indicates 
the mean value of [Ca/Fe] for Galactic GCs obtained from 
high-dispersion analyses (Gratton et al. 2004 and references 
therein). 



Our solutions for the a-element abundances of the 
VCC 1087 GCs are plotted in Figure U3J This is com- 
pared to the Galactic GC sample. Arguably, the most 
reliable "pure" a-capture element (i.e., that which is 
least affected by nucleosynthetic evolution and most eas- 
ily measured) is Ca (Gratton et al. 2004), and we plot the 
mean [Ca/Fe] ratio from high dispersion analyses in Fig- 
ure US For the Galactic GCs, ([E/Fe]) ~[Ca/Fe]=0.25, 
with dispersions of 0.16 and 0.11 dex respectively. The 
VCC 1087 GCs have ([E/Fe]) = 0.28 with a dispersion 
of 0.16 dex. Among these GCs there appears to be a 
weak correlation in the sense that [E/Fe] decreases with 
increasing [m/H], but this is does not have statistical 
significance. The position of VCC 1087 in Figure ITSlsug- 
gests a solar cv-abundance, which is consistent with that 
found by Geha et al. (2003). 



4. DISCUSSION 

With My = -17.8, VCC 1087 lies at the bright end 
of the dE luminosity function (e.g., Sandage, Binggeli & 
Tammann 1985). As such, specific conclusions for this 
one object are not necessarily generalizable to dEs as a 
class. This caveat notwithstanding, we discuss the find- 
ings of this study in the context of dE galaxy formation. 

Moore et al. (1998) and Mastropietro et al. (2005) pre- 
sented a series of numerical simulations which demon- 
strated that spiral galaxies can be transformed into 
spheroidal galaxies (i.e., dE/dSph types) in the cluster 
environment, "galaxy harassment". The combined ef- 
fects of the cluster tidal field and rapid encounters with 
cluster galaxies act to heat infalling progenitor spirals 
into, ultimately, prolate systems. Could VCC 1087 have 
been transformed in such a manner, and does its GC 
system offer any clues? 

One constraint on this picture is the number of GCs 
per unit (V-band) luminosity (S N ). Miller et al. (1998) 
measured Sn for 24 dE galaxies in the Virgo and Fornax 
clusters, and in the Leo group, 13 of which are of the 
dE,N class. These authors obtained a mean Sn of 6.5 ± 
1.2, and found a trend of decreasing Sn with increasing 
My. Miller et al. argued that luminous dE,N galaxies 
are unlikely to have evolved from faded irregular galaxies 
since locally irregulars exhibit Sn ~0.5-2. 5 We estimate 
a total GC population of 77±19 GCs in VCC 1087, giving 

5 Interestingly, Seth et al. (2004) have recently found that Virgo 
and Fornax cluster irregulars may have a larger spread in Sn than 
their field counterparts, with Sn ranging from ~ to 9. If con- 
firmed, this would allow for the transformation of cluster dlrrs to 
dEs (at least, that is, on the basis of an Sn argument). 
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Sn = 5.8 ±1.4 (adopting the slightly brighter value of 
My =-17.96 for VCC 1087 (Jerjen et al. 2004), we obtain 
Sn = 5.1 ± 1.2). These values are significantly higher 
than the typical Sn ~ 1 found for late- type spirals (e.g., 
Olsen et al. 2004; Chandar et al. 2004). In order to 
harbor ~ 80 GCs, a progenitor Sn = 1 spiral galaxy 
would need to have had My ~ —19.8. Any progenitor 
spiral must also have been essentially "pure disk" , since 
the densities of spiral bulges are too high to be strongly 
effected by harassment (Moore et al. 1998). 

The present-day Galaxy luminosity function in the 
Virgo cluster shows that, for My ~ —20, Irr galaxies 
are not represented, and Sd galaxies barely so (Binggeli 
et al. 1988). However, the peak of a Gaussian luminosity 
function for Sc galaxies does occur near this magnitude. 
An Sc galaxy which fell into the Virgo cluster ~ 5 Gyr 
ago (e.g., Conselice et al. 2001), and consequently lost 
its gas (via ram-pressure stripping; Gunn & Gott 1972), 
could plausibly fade by 1-2 mag in the V-band. Subse- 
quent mass loss through tidal encounters could fade the 
galaxy further, although this would presumably result in 
GCs being stripped also, lowering the Sn of the galaxy. 
Nevertheless, in terms of GC numbers, a faded Sc spiral 
could approach the S N of VCC 1087. 

By the same token, the high Sn of VCC 1087 is incon- 
sistent with the characteristic Sn ~ 2 observed for low- 
luminosity ellipticals (e.g., Ashman & Zepf 1998). Ad- 
mittedly, data is currently scarce on the Sx of such galax- 
ies, particularly regarding the question of the variation 
of Sn with environment. However, the present observa- 
tions sit uncomfortably with the notion that VCC 1087 
may be simply at the faint end of low-luminosity ellipti- 
cal galaxies (c.f., Miller et al. 1998). 

In principle, stellar kinematics provide constraints on 
dE formation models. Simulations suggest that galaxy 
harassment is more efficient at increasing velocity disper- 
sions (through bar formation and the transformation of 
circular to radial orbits) than disrupting rotation (Moore 
et al. 1998). Consequently, (v ro t / <Ji os )* should decrease 
due to the harassment process, but will remain non-zero 
except in cases where the progenitor spirals are strongly 
disrupted near to the cluster center (Mastropietro et 
al. 2005). Such severe disruption which may lead to non- 
rotating dEs is expected to be reflected in their observed 
morphologies and orbital characteristics. Observation- 
ally, Geha et al. (2003) found no obvious differences be- 
tween the stellar populations of a sample 13 non-rotating 
and 4 rotating dEs. Moreover, in the specific case of 
VCC 1087, Geha et al. (2003) found no evidence for ro- 
tation. 

This seems to present something of a challenge to the 
theoretical models. However, as discussed by Mastropi- 
etro et al. (2005), kinematical observations of dE galaxies 
are fairly sparse, and not without controversy. There are 
five galaxies in common between Geha et al. (2003) and 
the kinematical study of van Zee et al. (2004b). The 
latter study generally found higher maximum rotation 
velocities. Mastropietro et al. argue that this difference 
has arisen due to the fact that the van Zee et al. (2004b) 
observations covered a greater radial extent than those of 
Geha et al. (2003). Clearly further studies are required. 

In terms of GC kinematics the predictions are less clear 
since, to date, simulations have not explicitly taken them 
into account. Presumably, they will behave similarly to 



the dissipationless stellar component of the harassed spi- 
rals if they have the same density profile as the galaxy. 
Under the assumption that the azimuthal spatial distri- 
butions of GCs follows that of the galaxy stars, the GC 
systems of spheroids generally possess a significant de- 
gree of pressure support (Figure IBJ). However, the GC 
system of VCC 1087 appears to be rotationally flattened, 
which is possibly also the case for the only other GC sys- 
tem of a luminous dE studied to date, NGC 3115 DW1 
(Puzia et al. 2000). In contrast, the old halo Milky Way 
GCs, of which the MPCs in VCC 1087 are presumably 
analogs, show no significant rotation. We note that if in 
fact the VCC 1087 GCs represent the vestiges of a disk 
system, their rotation (i.e., an 'HI line-width' of ~ 210 km 
s _1 ) would place such a disk on the V-band Tully-Fisher re- 
lation at Mv ~— 18.7 (Sakai et al. 2000). If it is subsequently 
shown that the rotation in the GCs persists out to larger radii, 
this would correspondingly increase this magnitude estimate. 

The nucleus of VCC 1087 has a half-light radius (~ 3 pc) 
similar to the GCs, but appears ~ 0.1 mag redder. Tremaine, 
Ostriker & Spitzer (1975) first suggested that galaxy nuclei 
(in their case, that of M31) may be built from the merged 
remnants of massive GCs whose orbits decayed by dynamical 
friction. This idea becomes particularly appealing considering 
the evidence for a deficit of GCs, or a "core" , in the GC spatial 
density profiles in the inner regions of dE galaxies (Durrell 
et al. 1996a; Lotz et al. 2001; Section OHJ. Using Monte 
Carlo simulations, Lotz et al. (2001) found that the brightest 
nuclei of dE galaxies (Mv ~ —12) may well have been formed 
in such a manner, but the nuclei of fainter dEs were several 
magnitudes fainter than predicted (see also Miller et al. 1998). 
For VCC 1087, we obtain My ~ -11.4, which is a borderline 
case in the Lotz et al. (2001) simulations (see their figure 8). 

However, the color of the nucleus presents something of a 
problem for this scenario. As mentioned above, the nucleus 
has a color consistent with, or slightly redder than the pu- 
tative MRCs. This would suggest that the nucleus is solely 
constructed from decayed MRCs. This seems somewhat un- 
likely, considering that the VCC 1087 GC system consists 
predominantly of MPCs. It is possible that whatever lead to 
the formation of the MRCs formed the majority with a very 
high central concentration, leading to a shorter timescale for 
dynamical friction. Numerical simulations would be useful for 
investigating this issue. Interestingly, the Moore et al. (1998) 
simulations with gas also predict the existence of a possible 
nucleus. Gas sinks to the center of the galaxy in their simu- 
lations due to the combination of tidal torques and cooling, 
leading to a central density excess in their model galaxies. 

From the previous considerations, a plausible evolutionary 
scenario for VCC 1087 may have proceeded as follows: an 
Sc spiral with Mv ~ -20 and ~100 GCs fell into the Virgo 
cluster approximately 5 Gyr ago (the luminosity-weighted age 
of VCC 1087). Gas removal through ram- pressure truncated 
further star formation, and the galaxy stellar populations pas- 
sively evolved and faded by ~1.5 mag. Subsequent tidal inter- 
actions with cluster members and the cluster tidal field per- 
turbed and heated the galaxy, disrupting the disk and form- 
ing a low surface brightness spheroid. Possibly, some field 
stars and GCs where also stripped during this process. Any 
remaining gas in the system lost angular momentum from 
torques and cooled, leading to star formation in the central 
regions creating a solar metallicity nucleus, embedded in a 
diffuse envelope with a retinue of surviving GCs. 

5. SUMMARY 

We have performed a photometric and spectroscopic analy- 
sis of the GC system of the Virgo dE galaxy VCC 1087. This 
is the first such study of a dE galaxy in a cluster environ- 
ment. A total of 68 GCs are detected in HST/ ACS imaging 
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to (?475 ~ 26.5. Correcting for areal incompleteness and for 
the faint end of the luminosity function this gives a total GC 
population of 77 ± 19 GCs, corresponding to Sn = 5.8 ± 1.4 
in the V-band. The color-magnitude diagram of the GCs re- 
veals a pronounced blue peak ((/475 — Z850 ~ 0.90) and a tail 
of red GCs (g 4 75 - zsso ~ 1.17). Thus, VCC 1087 may have 
a metal-rich population of GCs similar to those seen in lu- 
minous ellipticals. Although poorly constrained due to small 
numbers, the spatial density profile of the GCs is surprisingly 
shallow, with a power-law exponent of a = —1.1 ± 0.4. This 
is similar to the spatial distribution of GCs seen in cD galax- 
ies, rather than those in low-luminosity ellipticals. The GC 
luminosity function (GCLF) appears log-normal, and has a 
turn-over magnitude very similar to those in luminous Virgo 
ellipticals (Strader et al. 2005). 

We investigated the kinematics and stellar populations of 
a subset of 12 GCs using Keck/LRIS spectroscopy. The GCs 
show significant rotation (V^ = 104 ± 35 km s _1 ). With a 
spatially-averaged velocity dispersion of trios = 29 ± 10 km 
s _1 , this suggests that this rotation is dynamically signifi- 
cant with (frot/cios)* > I- Comparing existing kinematical 
data for GC systems in early type galaxies suggests that the 
VCC 1087 GCs are among the most rotationally dominated. 
Considering both the Jeans equation and the tracer mass esti- 
mator (Evans et al. 2003) we derive V-band mass-to-light ra- 
tios (Ty) for the galaxy in the range 1 < Ty < 8 within ~ 6.5 
kpc of the galaxy center. These values are consistent with Ty 
of the galaxy stars, and can be understood purely in terms of 
a baryonic old, metal-rich stellar population. Stellar popula- 



tion analysis using Lick indices shows that the GCs are all old 
( J>10 Gyr), with a range of metallicities. We also estimate 
[a/Fe]~0.3, similar to the Galactic GC system. In contrast, 
the luminosity-weighted galaxy starlight has approximately 
solar metallicity, solar [cv/Fe], and is intermediate-aged (~ 4 
Gyr). 

We discuss the characteristics of VCC 1087 and its GC 
system in the context of the galaxy "harassment" scenario 
(e.g., Moore et al. 1998). We conclude that the observations 
are consistent with VCC 1087 being the remains of a faded, 
tidally perturbed Sc spiral, although this is probably not a 
unique interpretation of these data. 
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TABLE 1 

Basic Data for VCC 1087 Globular Clusters 



r^r 

ID 




o(J2UUU) 


9475 


Z850 


9475 — ^850 


77- 

Velocity 




(degrees) 


(degrees) 


(AB mag) 


(AB mag) 


(AB mag) 


(km s 1 ) 


GC1 


1 87 fl^fid-S 


J. J- . I oooo 


01 eoc; _l_ n nrj7 
il.OOO ZlZ U.KJKJI 


90 Qfifi 4- 0D7 
iu.yuu ziz u.uut 


n QiQ 4-n nm 

u.yiy zlz u.uuj- 


uuo zlz ±£ 


GC2 


187.07696 


11.77774 


21.950 ± 0.008 


21.141 ± 0.008 


0.809 ± 0.011 


783 ± 9 


GC3 


187.06023 


11.79187 


22.176 ± 0.011 


21.305 ± 0.014 


0.871 ± 0.018 


701 ± 22 


GC4 


187.05713 


11.77568 


22.209 ± 0.008 


21.408 ± 0.008 


0.801 ± 0.011 


695 ± 14 


GC5 


187.05296 


11.78424 


22.342 ± 0.009 


21.485 ± 0.009 


0.857 ± 0.012 


665 ± 10 


GC6 


187.06281 


11.79148 


22.397 ± 0.017 


21.550 ± 0.025 


0.847 ± 0.030 


647 ± 10 


GC9 


187.06194 


11.78850 


22.955 ± 0.035 


21.892 ± 0.044 


1.063 ± 0.056 


665 ± 35 


GC14 


187.08450 


11.80643 


22.892 ± 0.014 


22.008 ± 0.032 


0.884 ± 0.035 


745 ± 13 


GC16 


187.06472 


11.79208 


22.990 ± 0.016 


22.123 ± 0.017 


0.867 ± 0.023 


701 ± 34 


GC22 


187.06795 


11.78524 


23.440 ± 0.017 


22.359 ± 0.015 


1.081 ± 0.022 


631 ± 49 


GC25 


187.07155 


11.78644 


22.299 ± 0.009 


21.378 ± 0.008 


0.921 ± 0.012 


717 ± 31 


GC26 


187.07288 


11.78654 


23.713 ± 0.021 


22.584 ± 0.018 


1.129 ± 0.027 


698 ± 17 


Nucleus 


187.06185 


11.78996 


20.216 ± 0.003 


18.928 ± 0.003 


1.295 ± 0.004 


646 ± 30 a 


a Velocity from 


integrated light of nucleus+stellar envelope. 
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TABLE 2 
Kinematics of GC Systems 



ID 






(^/o"los)MPCs 




(f0/o"los)MRCs 




Ellipticity 






Source 




NGC 524 






0.75 




0.40 




0.10 




Beasley et al. 2004 




VCC 1087 






3.60 








0.26 




this work. 






NGC 1399 






0.07 




0.07 




0.00 




Richtler et al. 2004 




NGC 3115 






0.50 








0.60 




Kuntschner et al. 2002 




NGC 3115 DW1 




1.00 (1.36) a 








0.30 




Puzia et al. 2000 




NGC 4472 






0.27 




0.19 




0.19 




Cote et al. 2003; Zepf et al. 2000 


NGC 4486 






0.43 




0.44 




0.14 




Cote et al. 2001 




NGC 4594 






0.27 








0.46 




Bridges et al. 1997 




NGC 5128 






0.21 




0.30 




0.20 




Peng et al. 2004 




a Corrected for velocity uncertainties of 50 km s 1 and using the five 


inner cluster velocities 


given 


in Puzia ct al. (2000). 














TABLE 3 


















Lick/IDS indices and 


UNCERTAINTIES FOR VCC 


1087 GLOBULAR CLUSTERS. 






ID 




IT X 

tlop 


CNi 


CN 2 


Ca4227 


G4300 


H 7A 


Htf 


Fe4384 


Oa44oo 


Fe4531 




r X i 

[A] 


r X i 

[A] 


[mag] 


[mag] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


GC1 


4.07 


2.50 


-0.102 


-0.072 


0.28 


2.94 


0.05 


1.31 


0.32 


0.32 


1.96 


± 


0.24 


0.17 


0.007 


0.009 


0.14 


0.24 


0.25 


0.16 


0.37 


0.19 


0.29 


GC2 


2.94 


2.62 


-0.039 


-0.026 


0.46 


1.59 


-0.23 


1.55 


3.79 


0.14 


0.28 




0.21 


0.14 


0.006 


0.008 


0.11 


0.20 


0.22 


0.14 


0.31 


0.17 


0.25 


GC3 


1.61 


1.79 


-0.021 


0.005 


0.51 


3.49 


-0.18 


1.46 


1.86 


0.41 


1.25 




0.24 


0.16 


0.007 


0.009 


0.13 


0.24 


0.25 


0.16 


0.36 


0.19 


0.30 


GC4 


3.18 


2.63 


-0.089 


-0.063 


0.10 


0.85 


2.31 


2.57 


0.09 


0.25 


0.72 




0.25 


0.17 


0.007 


0.009 


0.13 


0.26 


0.25 


0.16 


0.40 


0.20 


0.29 


GC5 


2.82 


2.85 


-0.081 


-0.037 


-0.27 


2.14 


0.80 


1.52 


0.03 


0.92 


2.25 




0.28 


0.19 


0.009 


0.010 


0.16 


0.29 


0.28 


0.17 


0.44 


0.22 


0.32 


GC6 


2.25 


2.01 


-0.059 


-0.026 


0.38 


2.67 


-0.49 


1.23 


1.96 


0.59 


1.57 




0.32 


0.22 


0.009 


0.011 


0.17 


0.31 


0.32 


0.19 


0.45 


0.24 


0.38 


GC9 


-0.74 


0.87 


0.037 


0.027 


-0.58 


3.49 


-1.07 


0.23 


-0.76 


0.82 


0.87 




0.56 


0.31 


0.020 


0.023 


0.37 


0.67 


0.76 


0.50 


1.23 


0.66 


0.96 


GC14 


4.75 


2.60 


-0.103 


-0.066 


0.45 


2.89 


0.70 


1.70 


0.94 


0.71 


1.85 




0.39 


0.29 


0.012 


0.015 


0.23 


0.38 


0.40 


0.24 


0.59 


0.29 


0.45 


GC16 


2.92 


2.13 


-0.089 


-0.049 


0.03 


4.27 


-0.35 


1.60 


0.54 


-0.84 


0.65 




0.53 


0.38 


0.015 


0.019 


0.30 


0.52 


0.56 


0.34 


0.69 


0.43 


0.67 


GC22 


0.69 


0.86 


-0.105 


-0.116 


0.92 


-0.10 


0.08 


-0.18 


1.80 


0.66 


0.78 




0.66 


0.45 


0.019 


0.023 


0.29 


0.71 


0.63 


0.38 


0.92 


0.47 


0.69 


GC25 


-1.56 


1.57 


-0.024 


-0.011 


0.04 


7.06 


-6.64 


-0.50 


5.62 


1.61 


0.99 




0.92 


0.56 


0.027 


0.033 


0.49 


0.67 


0.93 


0.54 


1.13 


0.64 


1.01 


GC26 


0.53 


0.79 


-0.066 


-0.038 


-0.10 


3.89 


1.11 


2.39 


0.90 


0.40 


-0.51 




0.42 


0.29 


0.012 


0.015 


0.23 


0.41 


0.43 


0.26 


0.62 


0.33 


0.54 


Galaxy 


-1.28 


0.85 


-0.046 


-0.023 


1.69 


5.03 


-3.65 


-0.45 


5.24 


2.05 


3.25 




0.36 


0.24 


0.010 


0.011 


0.17 


0.31 


0.36 


0.22 


0.45 


0.23 


0.36 
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TABLE 4 

Lick/IDS indices and uncertainties for VCC 1087 globular clusters 



ID 


C 2 4668 


H/3 


Fe5015 


Mgi 


Mg 2 


Mg b 


Fc5270 


Fc5335 


Fe5406 




[A] 


[A] 


[A] 


[mag] 


[mag] 


[A] 


[A] 


[A] 


[A] 


GC1 


0.80 


2.07 


2.59 


0.049 


0.083 


0.83 


1.18 


0.97 


0.45 


± 


0.44 


0.16 


0.37 


0.004 


0.004 


0.18 


0.20 


0.23 


0.16 


GC2 


—2.81 


2.50 


1.92 


0.019 


0.057 


0.88 


1.23 


1.18 


0.24 




0.39 


0.14 


0.33 


0.003 


0.004 


0.17 


0.18 


0.20 


0.16 


GC3 


1.27 


2.11 


2.16 


0.037 


0.112 


1.37 


1.20 


1.10 


0.57 




0.41 


0.17 


0.35 


0.004 


0.005 


0.19 


0.19 


0.24 


0.16 


GC4 


-0.23 


2.79 


1.79 


0.011 


0.025 


0.07 


0.76 


0.79 


0.49 




0.48 


0.18 


0.36 


0.004 


0.005 


0.20 


0.23 


0.25 


0.19 


GC5 


-0.04 


2.49 


1.25 


0.030 


0.085 


0.59 


1.37 


0.88 


0.45 




0.48 


0.18 


0.40 


0.004 


0.005 


0.20 


0.22 


0.26 


0.20 


GC6 


0.79 


2.28 


2.42 


0.030 


0.096 


1.29 


1.79 


1.37 


1.15 




0.55 


0.21 


0.47 


0.005 


0.006 


0.23 


0.26 


0.32 


0.22 


GC9 


-0.86 


2.44 


0.57 


0.031 


0.134 


1.12 


1.43 


0.61 


-0.19 




1.33 


0.50 


1.05 


0.010 


0.011 


0.42 


0.49 


0.56 


0.43 


GC14 


-0.56 


2.86 


1.06 


0.024 


0.093 


1.62 


1.32 


1.89 


0.48 




U. I u 


U.ZD 


U.Oo 


u.uuo 


U.UU I 


n 9Q 


U.Ol 


U.OU 




GC16 


1.03 


1.84 


2.04 


0.031 


0.048 


0.16 


-0.21 


1.32 


-0.18 




1.01 


0.39 


0.86 


0.009 


0.010 


0.42 


0.47 


0.53 


0.41 


GC22 


-0.14 


1.47 


1.70 


0.027 


0.089 


1.44 


1.39 


2.86 


1.16 




1.08 


0.40 


0.85 


0.009 


0.011 


0.46 


0.50 


0.52 


0.39 


GC25 


0.49 


1.76 


3.61 


0.017 


0.099 


2.54 


1.41 


1.49 


1.27 




1.53 


0.58 


1.17 


0.014 


0.016 


0.63 


0.71 


0.80 


0.57 


GC26 


1.48 


2.11 


1.02 


0.020 


0.069 


0.41 


1.13 


1.05 


-0.04 




0.71 


0.32 


0.65 


0.007 


0.008 


0.34 


0.38 


0.42 


0.30 


Galaxy 


2.70 


2.07 


5.31 


0.069 


0.178 


2.82 


3.06 


2.10 


1.59 




0.55 


0.21 


0.45 


0.005 


0.006 


0.22 


0.24 


0.27 


0.20 



